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In this paper we consider some aspects of the relativistic dynamics of a cylindrical 
shell of counter rotating particles. In some sense these are the simplest systems 
with a physically acceptable matter content that display in a well defined sense 
an interaction with the radiative modes of the gravitational field. These systems 
have been analyzed previously, but in most cases resorting to approximations, or 
considering a particular form for the initial value data. Here we show that there exists 
a family of solutions where the space time inside the shell is flat and the equation of 
motion of the shell decouples completely from the gravitational modes. The motion 
of the shell is governed by an equation of the same form as that of a particle in a 
time independent one dimensional potential. We find that under appropriate initial 
conditions one can have collapsing, bounded periodic, and unbounded motions. We 
analyze and solve also the linearized equations that describe the dynamics of the 
system near a stable static solutions, keeping a regular interior. The surprising result 
here is that the motion of the shell is completely determined by the configuration 
of the radiative modes of the gravitational field. In particular, there are oscillating 
solutions for any chosen period, in contrast with the "approximately Newtonian 
plus small radiative corrections" motion expectation. We comment on the physical 
meaning of these results and provide some explicit examples. We also discuss the 
relation of our results to the initial value problem for the linearized dynamics of the 
shell. 

PACS numbers: 04.20.Jb,04.40.Dg 



I. INTRODUCTION 



In this paper we consider some aspects of the relativistic dynamics of a cylindrical shell 
of counter rotating particles. In some sense these are the simplest systems with a physically 
acceptable matter content that display in a well defined sense an interaction with the radia- 
tive modes of the gravitational field. The dynamics of these systems was analyzed originally 
by Apostolatos and Thorne but the evolution was considered in detail only over very 
short periods of time, and imposing aparticular form for the initial data, the "momentarily 
static radiation free" (MSRF) form [2], and the question of the general evolution in time 
of the system has remained largely unexplored. We notice that most of the literature that 
followed the work Apostolatos and Thorne has concentrated in the problem of collapse (see, 
for instance, [3[ and [H), and in general imposing particular forms for the fields, that may 
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include also some form of non gravitational radiation outside the shell (see, for instance, |5|, 
(fl or 0] ) . In a recent paper Hamity, Barraco and Cecere (sj , have considered again the rel- 
ativistic dynamics of these systems. In particular, since the system may have stable static 
configuration, and in the Newtonian limit small departures form the static configuration 
lead to periodic motions, it was expected that in the fully relativistic dynamics the inclusion 
of gravitational radiation modes should lead to a damping of these oscillations, through 
some form of "radiation reaction" . This expectation appears to be satisfied in the numerical 
solutions obtained in jsj. A closer analysis reveals, however, that the authors assumed an 
approximation where the back reaction of the radiative modes is essentially disregarded. 
This approximation would be justified if the coupling to the gravitational radiation modes 
had only a small effect on the dynamics of the shell. It turns out, however, as is shown in 
the present paper, that rather the opposite situation holds, and the dynamics is completely 
dominated by the behaviour of these modes. In fact we find that, in some sense, the coupling 
of the shell to the gravitational radiation modes is as strong as it can be, a remarkable fact 
that shows the dynamics of this system cannot be approximated by a Newtonian dynamics 
plus post - Newtonian corrections, as in the case of models where matter is confined to a 
bounded region. 

The plan of the paper is as follows. After setting up the problem in Section II, we show 
in Section III that there exists a family of solutions where the space time inside the shell 
is flat and the equation of motion of the shell decouples completely from the gravitational 
modes. The motion of the shell is governed by an equation of the same form as that of a 
particle in a time independent one dimensional potential. We find that under appropriate 
initial conditions one can have collapsing, bounded periodic, or unbounded motions. Next, 
in Section V we analyze the linearized equations that describe the dynamics of the system 
near a stable static solutions, keeping a regular interior. The surprising result here is that 
the motion of the shell is completely determined by the configuration of the radiative modes 
of the gravitational field. In particular, there are oscillating solutions for any chosen period, 
in contrast with the "approximately Newtonian plus small radiative corrections" motion 
expectation. Another interesting modes that appear here are the "anti resonances" discussed 
in Section VI. In Section VII we consider the general behaviuor of the periodic solutions, 
and in Section VIII their relation to the initial value problem for the linearized dynamics of 
the shell. We comment on the physical meaning of these results and provide some explicit 
examples. We also consider the role of the momentarily static and radiation free initial data 
of [H, in this context. Some closing comments are contained in Section IX. 



II. EQUATIONS OF MOTION 

We consider a spacetime M = M~ USU M + (M~ andM + are manifolds with boundary 
where the boundaries are identified with the 3-manifold S) with cylindrical symmetry where 
£ is the history of a hollow cylinder composed of counter-rotating particles of rest mass equal 
to unity; M~(M + ) is the vacuum interior (exterior) region of the cylinder. In the vacuum 
interior (M~) and exterior (M + ) of the shell, we introduce canonical cylindrical coordinates 
(t, r, z, 0). The metric takes the form 

ds 2 ± = e 2 ^" 2 ^ (dr 2 - dt 2 ± ) + e 2 ^dz 2 + e~ 2 ^r 2 d^ (1) 

Dropping the ± indices, the Einstein field equations in the empty space inside and outside 
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the shell are, 

1p,rr + -lff, r ~ ^ ,tt = (2) 

y >t = 2rij>,rtt > l,r = r[^, r ) 2 + ^,t) 2 } (3) 

We may interpret ip(r, t) as playing the role of a gravitational field whose static part is the 
analogue of the Newtonian potential. The time dependent solutions of (J2J) represent gravi- 
tational waves [||. Equation (J2J) is the integrability condition of Eqs. ([3]). The coordinates 
(z, 0, r) and the metric function ip are continuous across the shell E, while t and the metric 
function 7 are discontinuous. Smoothness of the spacetime geometry on the axis r = 
requires that 7 = and ip finite at r = 0. The junction conditions of M~ and M + through 
E require the continuity of the metric and specify the jump of the extrinsic curvature K ± 
compatible with the stress energy tensor on the shell. The induced metric on E is given by 

dsl = -dr 2 + e 2 ^dz 2 + e- 2 ^R 2 d<p 2 (4) 

Here V's(t) = i>+(R(i~), £+(t)) = i/)-(R(t), t_(r)). The evolution of the shell is character- 
ized by R(t), which is the radial coordinate r at the shell's location and r the proper time 
of an observer at rest on E. If we assume, as in [l|, and 0, that the shells is made up of 
equal mass counter rotating particles, the Einstein field equations on the shell may be put 
in the form, 

2A 

V£ - = 7=| : 2 . T2 (5) 

V R 2 + e 2 ^ s J 2 



4X\/R 2 + e 2 ^ J 2 
X + -X- = - VK +C J (6) 

R 

where the constants A and J are, respectively, the proper mass per unit Killing length of the 
cylinder and the angular momentum per unit mass of the particles. The other quantities in 
( )5|6|) are given by, 

dt± 



X ± = —± = + Ve- 2 fe-^) + R 2 (7) 
or 

fft = il$X ± + i/$R (8) 
where a dot indicates a r derivative, and we also have, 



d 2 R 
Ifr 2 



Rips - R 



R 2 ^; n X~ \R 2 X- J 2 e 2 ^X-X+ 

+ R 2 + e 2 ^J 2 ~ (R 2 + e 2 ^ J2)3/2 + ^2 + e 2^ J2) 

Equations (15TI6TI9]) , together with (|2|3|) determine the evolution of the shell and of the 
gravitational field to which it is coupled. The relevant functions: R(t), ip±{f-,t±), and 
7±(r, t±) appear satisfying a rather complex set of coupled ordinary and partial differential 
equations, with the boundary values for ip± and 7± at t = t±(r), r = R{r) directly coupled to 
the motion of the shell. Because of this complexity, the system was first analyzed in [I] onl 
to show some properties of the motion, although no solution was obtained, and later in p 
where, after introducing a second shell, mainly for technical reasons, an approximation that 
leads to an effective decoupling of (Q was used, to avoid considering the complex boundary 
problem that results for the wave equations (j2J) for ip±. Some full solutions of the problem 
are considered in the following Sections. 
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III. A RESTRICTED SET OF SOLUTIONS 

A full solution of the problem should provide the evolution of arbitrary initial data, 
satisfying the constraints imposed by the field equations. This is, clearly, a very complex 
problem. There is, however, a restriction on the set of solutions that while retaining its most 
interesting feature, namely, the coupling of the shell with radiative modes of the gravitational 
field, still simplifies considerably the system, allowing for a complete analysis of the resulting 
evolution and of its physical meaning. 



A. Static solutions 

We will first consider the static solutions for a shell of constant radius R, assuming an 
empty flat interior [H, (§]. In this case we may take 7_ = 0, ip_ = 0, implying tjj^ = 0, and 
for the exterior field we have, 

iP+(r) = -K\n(r/R) , 7 + (r) = 7o + k 2 \n(r/R) (10) 

then, 

X+(r) = e" 70 , X" = 1 (11) 
Since R = 0, tpY. = 0, and = 0, we find, 

J_ 2 

and, 



K = 2 h ' 7o = 21n[(i? 2 + 2J 2 )/it; 2 ] (12) 



A ~ (2J 2 + J R 2 ) 2 

This means that for any static solution we must have A < 0.15879..., (see 



B. Non static solutions with a flat interior 

The previous results indicate that, at least for the static case, we have solutions where the 
interior region of the shell is empty and flat. We notice that for a similar problem, namely 
a shell of counter rotating particles, but with spherical symmetry, we may have non static 
solutions where the radius of the shell changes in time, but the interior remains flat. In this 
case the spherical symmetry is crucial, as this implies that there are no radiative modes for 
the gravitational field. This is not the case for cylindrical symmetry, and, in general, one 
does not expect that in the non static case the interior will remain flat, because radiative 
gravitational modes, corresponding to a non static ip, will in general penetrate the interior 
region for, otherwise, the matching conditions would not be satisfied. Nevertheless, given 
the existence of the static solution with an empty flat interior, it is worthwhile to explore to 
what extent, if any, this condition can be generalized to a non static solution. We, therefore, 



assume again 7_ = 0, ij)- — 0, (implying ^ s = 0, ip^ n = 0, and X — y 1 + (R) 2 ), but place 
no restriction on either R(t) or and 7 + . The field equations are now (J2J), ([3]), and on 
the shell we have, 

2A 



X+ = \ 1 + {Rf - - (14) 
v ; VR 2 + J 2 
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Using this, and the fact that ip-£ = 0, we find, 



d 2 R _ (1 + (R) 2 )J 2 V 1 + (R) 2 (R 2 + 2J 2 ) 2 A 

dr 2 R{R 2 + J 2 ) R 2 (R 2 + J 2 f/ 2 ( ' 

The first surprising thing about this equation is that it contains no information on ift, 
and, therefore, it is an autonomous equation, completely decoupled from the gravitational 
mode. Equally unexpected is that it admits a simple first integral, given by, 



where C is a constant. This may also be written in the form 

(dR\ 2 [4A J 2 - XR 2 \n(R 2 + J 2 ) - 2CR 2 

\dr) + 4i? 2 ( J R 2 + J 2 ) 
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(17) 



and, therefore, the motion of the shell is identical to that of a particle of unit mass in the 
potential, 



V(R) = i 



[4 A J 2 - Ai? 2 ln(i? 2 + J 2 ) - 2C7? 2 ] 2 
4i? 2 ( J R 2 + J 2 ) 



(18) 



with vanishing total energy. Notice that C is not this energy and, therefore, the form of 
V(R) will be different for different solutions. Nevertheless, ( I17j) implies that if there are 
suitable choices of the parameters for which the potential V(r) has a negative minimum, 
the shell may execute a periodic motion. Let us first find the conditions under which this 
may happen. We look for equilibrium points (static solutions) where R = 0, and R = 0. 
Let R = R Q be that point, then, from ffT5l) . we have, 



_ j 2 r ^rJTj 2 

A - (R 2 + 2 J 2 ) 2 (19) 

which may also be considered as an equation for R , given A and J, 

To check for stable equilibrium points we set R(r) = R + £(r), replace in (fT5l) . and 
expand to first order in £(r). We find, 

^ [2^-(i? 2 + 2J 2 )J 2 ]J 2 ^ 

rfr 2 R 2 (R 2 + J 2 ) 2 (R 2 + 2J 2 f 1 J 



Then, the static solution will be stable for R 2 / J 2 > (1 + v / 17)/4, (i2o/ J > 1.1317...) and 
unstable otherwise (see also jsj). 

The somewhat complex form and dependence on its parameters of V(R) makes a general 
analysis of the possible motions based on (fT7|) rather difficult. We notice, however, that we 
have, 

V(R) = -(Xln{R) + C) 2 /2 + 0{R ) , R -)• oo (21) 

and 

V(R) = -2\ 2 R~ 2 + O(R ) , R -> (22) 
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FIG. 1: Plots of V(r) as a function of ln(r). We have taken J = 1, and A = 0.1 in all cases. The 
solid curve, for C = —1.3, contains a region where periodic motions are possible. This region is 
absent for both the dotted [C = —1.1), and the dashed (C = —1.5) curves. In all cases we have a 
collapsing region, for sufficiently small r, and unbounded region, for sufficiently large r. 



and therefore we have unbounded motions for sufficiently large R and collapsing motions 
for sufficiently small R. Moreover, the equation, 

4 A J 2 - XR 2 \n(R 2 + J 2 ) - 2CR 2 = (23) 

has a real root with R = R m > for any real J and C, and A > 0. But for R = R m 
we have V(R) = 1/2, which is also the maximum possible value of V(R), and, therefore, 
the collapsing and unbounded motion regions are separated at least by a "forbidden" gap. 
Depending on the values of the parameters, V(r) may contain two "forbidden" gaps, where 
V(r) > 0, and periodic motions are possible in the region between these gaps. Figure 1 
provides some explicit examples of these cases. They will be explored in more detail in the 
next Sections. 

To close this Subsection we remark also that the evolution equation ()15p has the following 
scaling property: if we introduce the function R(rj), such that, 

R(rj) = -R(Jrj) (24) 



J 



we have 



dv 2 (& + i y /2 R2 r(r 2 + i) 

and, therefore, all the types of motions, up to scalings, are determined by the (adimensional) 
parameter A. 
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C. Compatibility with the field equations 

So far we have only considered Eq. ()15p . The full set of field equations includes also ip, 
7 and the junction conditions, and there is, a priori, no guarantee that the only solutions of 
( TT5l) compatible with these are the static ones. In particular, the condition = implies, 

^(R(r),t + (r))R + ^(R(r),t + (r))X + = (26) 

which, together with ([8]), and (TL4")) . determine ^+(i?(r), t + (r)) and ^(R(r), t + (r)) in terms 
of R(t). After some simplifications, and using ( TT7j) . we find, 



#(i2(r),t + (r)) 



A[4AJ 2 + 8Ai? 2 + \R 2 \n(R 2 + J 2 ) + 2CR 2 } 
R(R 2 + J 2 ) [1 + 4A 2 ln(i? 2 + J 2 ) + 8AC + 16A 2 ] 



t/j + t (R(T),t + (r)) = - — (27) 

'* V V ' V " VR 2 + J 2 [1+ 4A 2 \n(R 2 + J 2 ) + 8AC + 16A 2 ] V 7 



Similarly, from (1141) . we have 



7 + (i?(r), t + (r)) = ln[l + 4A 2 ln(R 2 + J 2 ) + 8AC + 16A 2 ] (28) 



From this equation we may compute, 

dry* 



dr 



1+R + 7t^ + (29) 



and we can check that if we replace ([3]) on the right, and then use ([27)1 . we get the same 
expression as that obtained by computing the left hand side of (129]) using ([28~]) . We conclude 
that the restriction to a flat interior is compatible with the dynamics of 7 on the shell, even 
in the non stationary case. 

Now we could compute in principle ip (and then 7) outside the shell. We notice that, 
provided R(t) satisfies some suitable conditions, to be considered below, if we only imposed 
ipn = 0, since R{t) is given, then we would get for ip a wave equation with a well defined 
boundary condition . But in this case both ifj^ and ip^ are given in the boundary, and it is 
not clear that in this case we may get any non trivial solution. To analyze this problem we 
notice that in ()2]) we may consider r as the "time" variable, and t as the "space" variable 
(Io| . as shown in Figure 2. Then, the problem can be posed as that of finding the evolution 
(in r) of ip, for initial data (V's = 0, ip >r ) on the (one dimensional) surface S, 

t = v , r = R{r{rj)) (30) 

where rj is a parameter, and r(t) is obtained by inverting (I14p . The problem is well posed 
provided that S is a Cauchy surface, and this requires that the tangent vector to S be "space 
like", that is, (dt/dr]) 2 - (dr/dr/) 2 > 0, or 

(X+) 2 >(7?) 2 (31) 

But, from (J7J), this is always satisfied. We might, therefore, conclude that the field equations 
have solutions for any R(t) that is a solution of (115)) . However, we must also require that 
the data ip tr (R(r)) be non singular, but, as can be seen from ( 127)) . this may not always be 
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FIG. 2: Evolving ip + {t,r) from "initial data" = 0,?/v) on the curve (r = R(r),t = t + (r)), 
represented by the thick curve in the Figure. tp + (t,r) is evolved into the region above this curve. 



the case, because the factor [1 + 4A 2 ln(i? 2 + J 2 ) + 8AC + 16A 2 ] in the denominator in (|2 
may vanish for some finite R ^ 0. We notice here that solving (J23l) for In^R^ + J 2 ) we get 



1 + 4A 2 m(i4 + J 2 ) + 8AC + 16A 2 = 1 + 16 



{Rj + J 2 ) A 2 

p 2 



(32) 



Since the left hand side is a monotonically increasing function of R, this implies that the 
denominator is always positive for R > R m , and, therefore, this problem does not arise for 
the unbounded solutions of the previous Subsection. Periodic motions are only possible if 
the potential has, besides that for R = R m with V(R m ) = 1/2, another maximum for say 
R = R p , with V(R P ) > 0. For this maximum we would have, 



dV 

Jr 



FiF 2 



R — Rr> 



where, 



AR 3 (R 2 + J 2 ) 2 
F 1 = XR 2 \n{R 2 + J 2 ) + ICR 2 - 4 A J 2 







(33) 



(34) 



and, 



F 2 = R 2 X In (R 2 + J 2 ) J 2 + 2 A (2 J 4 + R A ) + 2 R 2 J 2 (C + 4 A) (35) 
The first factor in (1331 . Fi, vanishes for i? = i? m . We can check that this is its only zero 



for R > by noticing that, 



dF l 



dR 



R—Rn 



X (R m 2 + 2J 2 Y 

Rm (Rm 2 + J 2 ) 



> 



(36) 



and, therefore, the equation i*\ = can have only one root for R > 0. Then, any extremum 
of V^(i?) other than that for R = R m must come from the vanishing of the second factor, F 2 . 
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It can be shown that F 2 can have zero, one or two roots depending on the parameters and 
that each root R p must satisfy R p < R m . If F 2 has no root or only one, V(r) has only one 
maximum (at R m , the root of F 2 being a saddle point) and there are no periodic motions. 
If F 2 has two roots, V(r) has two maxima and a minimum between them (a maximum and 
a minimum at the roots of F 2 and another maximum at Rm). In this case we can have 
periodic motions only if V > at the first maximum and V < at the minimum, and this 
depends non trivially on the choice of parameters. We can, nevertheless, obtain a useful 
result as follows. At a root R p of F 2 we have, 

In (R p 2 + J 2 ) = ~2 2j4 + f/ j2 ARp2j2 ~ 2y (37) 



and, therefore, we may write, 



Then, if V(R P ) is a positive maximum we must have, 

a2 = rj J*R P 2 (R p 2 + J 2 ) 
(i? p 2 + 2J2) 4 

with < rj < 1. Replacing in the problematic factor in the denominator of (1271) we have, 

,„ , , _ A x9 J 2 (i?„ 2 + J 2 ) (iL 4 + 2R P 2 J 2 + 2J 4 ) ?7 / 

1 + 4A 2 In (iL 2 + J ) + 8CX + 16A 2 = 1 - 8 — K — J -L* \ '— 40 

(R P 2 + 2J 2 ) 

The right hand side of (140]) is a linear function of 77. It is equal to 1 for 77 = and to 
R p 8 / (R p 2 + J 2 ) 4 for 77 = 1. Since the left hand side of (1401) is a monotonic function of i?, it is 
positive at any positive maximum of V(R) and the oscillating solutions always take place at 
R > R p (after the first positive maximum), we conclude that (1271) is well defined and finite 
for any periodic motion. 

As indicated, there are also collapsing solutions with a flat interior. In this case we find 
that ip + t (R(t), t + (r)) is singular either at some finite R or at R = 0. It can be checked 
that these singularities occur for finite r and t^. Since the solutions are symmetric in r, 
this implies that the evolution has singularities both at some finite time in the past and in 
the future, and, therefore, by causality they extend only to some bounded region in r. We 
do not analyze further these solutions as they do not seem to be physically interesting. 

In the following Sections we consider linearized solutions corresponding to infinitesimally 
small departures from the static stable solutions, both for flat and for empty regular interiors. 



IV. LINEARIZED PERIODIC SOLUTIONS WITH A FLAT INTERIOR 

Let us assume that R = Rq corresponds, for some suitable J, to a stable static solution 
with A given by (ITU]) . For this solution we have ip- — 0, ips — 0, and 7_ = 0, with ip + and 
7 + given by (I10|I12I) . We consider now a perturbation of the static solution such that the 
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interior remains flat. This means that we keep = 0, ip-z — 0, and 7_ = 0, but for the 
other dynamic variables we introduce now a time dependence by setting, 



R(r) = Ro + ^r) 
ip + (r,t+) = -Kln(r/Ro) +ip(r,t + ) 
7+(r,t+) = 7 + K 2 \n(r/R ) + 7(r,t + ) 



(41) 



with k and 70 given by (fT2l) with R = R , and consider the linearized field equations that 
result from expanding to first order in £, if; and 7. To this order £ satisfies (120]) . If we define, 



[2R 4 - (R 2 + 2J 2 )J 2 ]J 2 



Q 2 
"0 



(42) 



R 2 (R 2 + Jiy(R 2 + 2Ji) 
the solution of (1201) can be written as, 

£(r) = £oe inor (43) 
In accordance with (EJ), and (JHj), and expanding to first order in £ > we have, 



4J%e in ° T 
dr ~ {R 2 + 2 J 2 ) 2 ' (# 2 + 2 J 2 ) 2 i? ( 



+ 



(44) 



Actually, the last term on the right in (|44p contributes in all relevant equations only to 
second order, and, therefore, we may set, 



Rq 



(R 2 + 2 J 2 ) 2 



(45) 



when appropriate. Similarly, we may set £_ = r. We also define, for convenience, 

(R 2 + 2J 2 ) 2 ^ 
R 



(46) 



We look now for solutions of ijj and 7 with the same periodicity as R(t). On account of 
and (145]) the general solution for ip will be then of the form, 



2J 2 
"ft 2 



In 



r 

R 



+ (A 2 J (n 2 r) + B 2 Y (n 2 r)) e iQ2t + 



(47) 



where A 2 and B 2 are constants, and Jo, and Y are Bessel functions. Then the junction 
conditions on the shell, and the condition ips = are satisfied (to first order in £ ) if> 



A 2 = -txJ 
Bo = irJ 2 



J 2 (AJ 2 R 2 + AJ 4 -R 4 
(R 2 + .P)Rl 



-Y^Ro) 



Q {R 2 + 2J 



2\2 



R% 



£0 



t2\2 



(fig + 

Similarly, again to first order in £ 0) we find 



7+(^+> r ) = 7o + 



4J 4 



Kq 



4J 2 
"ff 2 " 



(A 2 Jo(f2 2 r-) + 5 2 y (n 2 r))e 



(48) 



(49) 
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where, 

7o = ,JM±^!l (50) 

Summarizing, we see that given appropiate values of Rq and J, we can find a complete 
solution, at the linearized level, where both the motion of the shell and the radiative modes 
of the fields are periodic in their respective times. For this type of solutions the period is 
a definite function of i? an d J, in correspondence with the idea of a "perturbation" of a 
stable equilibrium static configuration, characterized by Rq and J, with the departure from 
equilibrium being given by the arbitrarily small parameter £ - Finally we remark that in 
the limit J 2 « Rq we have, 

2 J 2 

that is, Qq approaches the value corresponding to small oscillations of the shell in the New- 
tonian limit. At first sight it would appear that this should be the natural frequency of 
oscillation of the shell, and that the effect of the coupling to the gravitational radiation 
modes should introduce only a small departure, such as damping, from the Newtonian case. 
However, as we shall show in the next section, this is only a special case resulting from the 
assumption of a flat interior, and the behaviour of the system is in general quite different 
from this expectation. 



V. LINEARIZED PERIODIC SOLUTIONS WITH A REGULAR INTERIOR 

We consider now the more general situation where the interior region is empty but may 
contain gravitational radiation, imposing only the condition of regularity on the symmetry 
axis r = 0. We then set, 

V>-(<-,0 = A 1 J (fi_r)e iQ -'- (52) 
Restricting again to linearized order we may set, 

T _(t_,r) = (53) 

and, therefore, also to the appropriate order, we may also set, 

*_(t)=t (54) 

We assume again a perturbation around a stable equilibrium configuration characterized by 
Ro and J. We therefore take, 

R(T) = R + ^ nr (55) 

and, 

MU, r) = -j^- In (- + (A 2 J (Q 2 r) + B 2 Y (Q 2 r)) (56) 

where A 2 and B 2 are constants, considered to be of first order. To this order we then have, 

4/4 / r \ 4/2 

7+(f+,r) = 70 + ^hi — - (A 2 J (Q 2 r) + B 2 Y (Q 2 r)) e^ 4 + (57) 
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where 70 is given by (1501) . A long calculation then shows that consistency at first order of 
the equations requires = Q, 



and, 



(R 2 + 2 J 



2\2 



-fi = e 7o ft 



Rq 



-T 



e -7o r 



(58) 



(59) 



{Rl + 2 J 2 ) 2 

Replacing now in (jSJ), (El), and ([9]), and expanding to first order, we find a set of three linear 
independent equations for Ai, A2, B 2 , and £o- It turns out that a convenient way of handling 
this system is to introduce a new parameter a by the definition, 



R 2 (R 2 + 2J 2 ){R 2 + J 2 ) 2 Q 2 + J 2 (2J 4 - 2i? 4 + R2J 2 ) 
A\ — — „ n — : a 



q 2 r 2 



R 2 (R 2 + 2J 2 )(R 2 + J 2 ) 2 (n 2 -tt 2 ) 



n 2 R 2 



where Qq is given by (|42p . We then have, 



(go(g + 2j 2 )(r 2 + j^nj^nRp) - j 2 (2j 2 + 3g|)Jo(ago)) gj 

fi 2 i? 2 + 1 



and, 



(60) 



(61) 



A, 



h (i? 2 + 2 J 2 ) 2 (V + J 2 ) V - J 2 (J 2 + 2i*g)) Fx (fi 2J R ) 
+2J 2 F (V 2 R ) #0 (^0 + j2 ) (( 4j4 + 6 ^o j2 + R o) V 2 - 2 J 2 )] J (ORo) 
- V (fi 2 i? ) ((^ + 2J 2 ) 2 (i? 2 (i? 2 + J 2 ) 2 VL 2 - J 4 ) - 2J 2 ^) 

tt(R 2 + 2J 2 )a 



-2J 2 R n (R 2 + J 2 ) (i? 2 + 2 J 2 ) 2 Y 1 (n 2 R ) QJi (SlRo) 



2Rl{R 2 Q? + 1) 



(62) 



B 2 = - [(-2 J 2 Y (n 2 R ) R (Ro 2 + J 2 ) ( (4 J 4 + 6 i? V 2 + #o 4 ) ^ 2 - 2 J 2 ) 

-fi (i? 2 + 2 J 2 ) 2 (i? 2 (i?o 2 + ^ 2 ) 2 ^ 2 - J 2 (J 2 + 2 .Ro 2 )) *i (O 2 i?o)) Jo (fii2o) 
+Jx (Oi? ) (*o (fi 2 i? ) ((V + 2 J 2 ) 2 (r 2 (V + J 2 ) 2 « 2 - J 4 ) - 2 J 2 R { 



-2 J 2 Ron (Ro 2 + J 2 ) (Ro 2 + 2 J 2 ) 2 Y x (n 2 R ))] ^^^TO ( 63 ) 

The main reason for displaying these, at first sight, not very illuminating expressions for 
Ai, £0, A 2 and B 2 is that they explicitly show that given R and J corresponding to some 
stable equilibrium configuration, i.e., to some real value for Qq, we have non trivial periodic 
solutions for the linearized perturbations for every value of Q. Thus, we reach the unexpected 
result that, at least perturbatively, we cannot ascribe a particular period to motions close to 
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the stationary solution, as happens in the corresponding Newtonian dynamics. The period 
of the motion can be arbitrary, depending entirely on the field configuration. From a more 
physical point of view, this can be interpreted by noticing that as the radius of the shell 
changes, the change in the static part of the field (the ln(r) terms) is of the same order of 
magnitude as the radiating part of the field that this motion generates. Thus, as the shell 
moves away from its stationary configuration, the motion is driven to essentially similar 
extents by the static and the dynamic parts of the gravitational field. In some sense then, 
the coupling of the shell to the gravitational radiation modes is as strong as it can be, a 
remarkable fact that shows that the dynamics of this system cannot be approximated by a 
Newtonian dynamics plus post - Newtonian corrections, as in the case of some more realistic 
models, where matter is confined to a bounded region. 

We have already analyzed the special case where the field inside the shell vanishes, and 
found that this is possible only for a particular value of Q, which, in the context of this 
more general analysis, corresponds to the particular solution where A\ = 0. In fact it is 
straightforward to show that the solution for Q = Q Q reduces precisely to that of the previous 
Section. But there is also, for instance, a particular set of solutions that display a different 
type of unexpected behaviour. We may call these "anti resonances" . They are described in 
the next section. 



VI. ANTI-RESONANCES 

As the shell evolves in time, its physical radius is given by Ro(t) exp(— ^s( r ))- F° r 
perturbations around an equilibrium point, to linear order we then have, 

Ro(r)e-^ =R + (& - A 1 R J (QR )) e lQT (64) 

If we use now ([6"0]) and (l6Tj) . we get, 

^-A 1 R J (QR ) = [R 3 nj 1 (nR )-(j 2 + R 2 n 2 (R 2 + J 2 ))Jo(nR )] 

R (Ro 2 + 2J 2 ) (R 2 + J 2 )a 
x i '- — (65) 

R 2 n 2 + i v ; 

This implies that the physical radius of the shell remains constant (to first order), and hence 
we have an anti-resonance, if Q is a solution of the equation, 

R 3 n Ji (ft Rq) ~ {J 2 + Ro 2 V< 2 (R 2 + J 2 ) ) Jo (ft Ro) = (66) 



It is easy to check that (166|) has an infinite sequence of solutions. Again it is remarkable that 
for these frequencies the effects of the inner an outer radiation modes exactly compensate 
each other and the shell remains motionless. 



VII. THE GENERAL BEHAVIOUR OF THE PERIODIC SOLUTIONS 

The linearized solutions found in the previous sections have in common the desirable 
feature that they contain a parameter that can be made arbitrarily small, and thus they 
approach arbitrarily closely the static solution. At least this is true for finite values of r. In 
fact, looking at the form of (}l"Tf) and (|4"5j) we see that for large r the solution appears to be 
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dominated by the static ln(r) terms, and, therefore, that the solutions approach the static 
unperturbed background for large r. A more accurate geometrical picture can be obtained 
by considering, e. g., the Kretschmann invariant K, for large r. Using the forms (l47|l and 
we find, 



A" 



16k(1 + k)(1 + 2k) 

e 70_R- 4K ( 1+K ) r 4+4 K (l+ K ) 



cos(f2!t H 



2fi 3 r 3 



7T 



X 



A\ sin {^l\r + — J — A 2 cos (flir + ^ 



O 



,1/2 



r 4+4K(l+re) 



(67) 



while for the background static metric we have, 



K 



16k 2 (1 



« 2 )(1 



e 47o_^- 4K ( 1+K ) r 4+4K(l+K) 



(68) 



Thus, although in both cases K — > for large r, in the perturbed case K is a factor of order 
r 3 / 2 larger than in the static case, so this appears to indicate a larger and larger departure 
between the perturbed and perturbed solutions as r — > oo. The consequences and meaning 
of this departure are not clear. For instance, in the flat interior case, where we have the 
same behaviour for K, we have shown that there are non perturbative periodic solutions, 
and the linearized solutions should approach those, and therefore, at least in that sense, the 
behaviour f)67p would be compatible with a perturbative treatment. This, however, is not 
entirely correct. The reason is that if we attempt to solve the field equations for ip + and 7+ 
to second order in the periodic terms, 7+ acquires terms of order r (rather than order r° as 
in the first order terms), and, if these are included in K the difference between the solutions 
is now of order r 3 . 

There is nevertheless, another way to look at these solutions. As indicated, they are 
indeed close to the static solution provided r is not too large. We notice that the equations 
for ip + and 7+ are local and causal. In particular, d r y + /dt + vanishes if dip + /dt + = 0. We 
may, therefore, consider the solution up to some large value of r, say r b » R , cut off 
the periodic part for r > r^, and use this configuration as initial data for the system. On 
account of causality, the shell will then oscillate periodically for a time of the order of 
rb, while preserving the asymptotic structure of the static solution, so that, in principle, we 
could have solutions that are periodic for a time that is long as compared with the oscillation 
period. 

There is still another way to look at the linearized solutions that is explored in the next 
section. 



VIII. THE INITIAL VALUE PROBLEM 

An important problem related to the system under discussion is the following. Suppose 
we have initial data that differs slightly from that corresponding to the static solution. We 
then expect that the evolution of that data will remain close to static solution and that, 
therefore, a linearized treatment should be adequate. We notice at this point that a linear 
superposition of linearized periodic solutions will also be a linearized solution, although no 
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longer periodic. In fact, we may generalize this idea and write, 

roo 

ip+{t + ,r) = -n\n(r/R ) + / [a 2 (fi+) J (fi+r) + b 2 {n + )Y {n+r)} e in + t +d^t+ 

Jo 

roo 

-y+(t+,r) = 7o + ^ 2 ln(r/ J R )-2K / [a 2 (fi+) J (fi+r) + b 2 (n + )Y (n+r)] e* n+t+ dn + 

Jo 

roo 

tjj-(t_,r) = / ai(fi_)J (fi_r)e <n -*-dfi_ (69) 
Jo 

7_ = 

/■oo 

i?( r ) = R + t(n)e iQr dn 

Jo 

where the coefficients a\, a 2 , b 2 , and £ are complex functions of their arguments, and, as 
usual, it is understood that we take the real part of the right hand side of ( 169]) . Since each 
value of Q is independent of the others, we may use the results of the previous Section to 
solve for a 2 , b 2 , £ and a\ in terms of a different a, for each Q, and it is clear that we may 
choose a to be an arbitrary complex function Q. In particular, considering the asymptotic 
behaviour for large r of the Bessel functions Jo and Yq, we see that with an appropriate fall 
off for a(Q) as Q — > oo we may control the corresponding fall off for large r of and 7+, 
because the dependence on Q of both a 2 and b 2 is related by linearity to that of a. Since 
the expressions for and 7+ have the form of Fourier - Bessel transforms, in this case 
the radiative parts should also fall off for large \t\, and these expressions might represent 
a situation where for large negative t the shell is stationary, being subsequently perturbed 
by an incoming gravitational radiation pulse, which eventually rebounds, leaving the shell 
again in a stationary state. 

Although the above reasoning is correct, it is not clear how we can use it to solve the 
initial value problem for our system. To begin with, assuming that, e.g., ip + (0,r) is given, 
since the range of r is not < r < 00, and we cannot impose a priori boundary conditions 
for r = R , there appears to be no well defined procedure for inverting (|69|) and computing, 
say, a 2 and b 2 . Nevertheless, since a, and therefore, a 2 and b 2 , are complex, we actually 
have two arbitrary real functions of Q + at our disposal for the construction of V'+(0, r), and, 
therefore, make the system satisfy arbitrary initial data. We notice, however, that once a(Q) 
is given, not only a 2 and b 2 are fixed, but also a\, and therefore, the data inside the shell, 
which should, from causality, be independent of that outside the shell. The answer to this 
conundrum is that the expression for ifj + (Q,r) on the right in (1691) is overcomplete, because 
we only require if) + (Q,r) in the range _R < r < 00, and that leaves the range < r < Rq 
arbitrary. Since the expression in fl69|) actually defines ^+(0, r) also in the region < r < Ro, 
there must be an infinite set of functions a 2 and b 2 that reproduce the data in Rq < r < 00, 
so in principle there is room for arbitrary data ij)-(0, r) in < r < Rq. Again, although this 
seems plausible, we do not have a proof of its validity. The difficulty here is the lack of a self 
adjoint formulation for the initial value formulation of the moving boundary problem posed 



by the dynamics of our system. This problem will be considered in detail elsewhere [11 . 

To illustrate the points considered in this Section, we include as an example, the case of 
an incoming pulse, its interaction with the shell, and eventual rebound after this interaction. 
In this example we set Ro — 4, J — 1, which implies A = 0.0509..., and VLq = 0.0770.... We 
also set, 

, R 3 (n 2 R Q 2 + l)Qe-W- 2 ') 2 
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FIG. 3: Plot of £(t) = R(t) — Ro in the region —10 < r < 10. The region r < is dominated by 
the incoming pulse, while that for r > is dominated by the outgoing pulse, resulting from the 
rebound of the pulse on the symmetry axis. The shell is essentially in its equilibrium radius Ro for 
either r << —15 and r >> 15. 

with Q = 10 -5 . Replacing in (1691) we obtain explicit expressions for the dynamic variables 
of the problem, from which we can view the evolution of the system. Details are given in 
Figs. 3, 4 and 5. In Figure 3 we have a plot of R(t) — Ro as a function of r, showing the 
incoming pulse region, for r < 0, and the outgoing pulse region for r > 0. The shell is 
essentially in its equilibrium radius Ro for either r << —10 and r » 10. Figure 4 is a plot 
of ijj_(t,r) in the region < r < R = 4, —10 < t < 10. We notice the propagation of the 
incoming pulse, an intermediate interference zone, formed by incoming and outgoing waves, 
and the eventual fall off of the pulse as it propagates outside. In Figure 5 we have a plot 
of i/j + (t,r) in the region R Q = 4 < r < 15, —10 < t < 10. We can see the propagation of 
the incoming pulse towards the shell, a zone near t = where most of the pulse has gone 
through the shell, and its rebound and propagation away from the shell, for t > 0. 

In closing this Section we may ask what is the relation of this construction, and its implied 
initial data, to the "momentarily static and radiation free" initial data of [lj]. (See also [l2[ 
for a different analysis of the meaning of this type of data). To understand this we consider 
again the general equations. Without loss of generality we may choose r, t_, and t + such 
that t = corresponds to t + = t_ = 0, and then, for the "momentarily static and radiation 
free" initial data on the surfaces t + = £_ = r = 0, we would have, 



7_(t_,r) 



R(r) 



R 1 + R 2 r 2 /2 + . . . 
Pi{r) t 2 _ + ... 
Gi(r) 

Q 1 Hr/R 1 ) + P 2 (r) t 2 + + ... 
Q 2 1 ln(r/R 1 ) + Q 2 + G 2 (r) t\ + ... 



(71) 
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FIG. 4: Plot of ip_(t,r) in the region < r < R = 4, —10 < t < 10. We notice the propagation 
of the incoming pulse, an intermediate interference zone, formed by incoming and outgoing waves, 
and the eventual fall off of the pulse as it propagates outside. 
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FIG. 5: Plot of m the region 4 < r < 15, —20 < t < 20. We can see the propagation of 

the incoming pulse towards the shell, a zone near t = where most of the pulse has gone through 
the shell, and the rebound and propagation of the pulse away from the shell, for t > 0. 
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where Pi, P 2) Gi, and G 2 are some functions of r that are determined by the evolution, 
and whose explicit form is not relevant, and dots indicate higher order (in r and t±) terms. 
Replacing these expressions in the general equations of motion we find that we must set, 



Qi 



2\Rx' 



4 A [Ri 2 + J 2 ) - Riy/R^ + J 2 



o, = -^v^ ±g, (72) 



R2 



Ri 

(2 J 2 + i?x 2 ) 2 A - J 2 y/R 1 2 + J 2 R 1 
2 I T9\3/2 D 2 



(i?! 2 + J 2 )^^ 2 

We recall now that for a static solution we have, 

= JW + J 2 i?Q 
(2J 2 + i? 2 ) 2 

where i?o is the equilibrium radius. For small departures from this radius we may set, 

Ri = R + Z (74) 

where £ is constant. Replacing ( 1731) and (1741) in ( 1721) . and expanding to first order in £ we 
find, 

J 2 J 4 (-R 4 + 4 J 2 /^ 2 + 4 J 4 ) 
g 2 = -4 In (R ) + 2 In (2 J 2 + i? 2 ) - 4 ^£ + O (£ 2 ) (75) 



[2J 2 + J R 2 ) (Ro 2 + J 2 ) R 



o 2 



Thus, for sufficiently small departures from the static stable solution the "momentarily 
static and radiation free" initial data corresponds to adding a constant to iJ)+(0, r), while 
leaving ^_(0,r) = 0. As indicated already, it appears possible, at least in principle, to find 
coefficients in (|69|) that correspond to this initial data, and that could be used to study its 
evolution, but this problem has not yet been solved. As a final comment, we notice that 
this does not correspond to perturbations of essentially compact support, and therefore, one 
would expect that the corresponding integrals would contain some singular expression. This 
may not be a problem, because to analyze the motion for finite times near this initial data 
we may simply cut off the perturbation for large r, which, from causality, would not modify 
the solution for times of the order of r. We, nevertheless, refer to [l2[ for a different analysis 
of this problem. 



IX. COMMENTS 



In this paper we have presented an analysis of the dynamics of a self gravitating cylin- 
drical thin shell of counter rotating dust particles. This analysis provides several new and 
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to a certain extent unexpected results. In particular we show that there exists a family 
of solutions where the interior of the shell remains flat at all times. For this family the 
equation of motion for the radius of the shell decouples from the radiative modes. We find 
a first integral for this equation and show it to be equivalent to that of a particle in a one 
dimensional time independent potential for a certain value of its total energy. Depending 
on the constants of the motion we have collapsing, periodic or unbounded solutions. We 
further analyze under what conditions these solutions are consistent with the field equations 
for the gravitational modes, and show that there are consistent periodic solutions for the 
full system. We consider next the dynamics of the system close to a stable static solution 
in the linearized approximation, where we assume that the interior is regular. The first 
unexpected result is that we have non trivial solutions for any possible frequency, and that 
all these modes are stable. The only role played by the Newtonian frequency (corresponding 
to Newtonian dynamics of the shell) is that it is the only frequency for which the interior 
is flat. We thus reach the conclusion that the system has no "natural" oscillating frequency 
that would be slightly modified by the coupling to the radiative modes, but, rather, we have 
a system where this coupling is "as strong as it can be" , fully determining the behaviour of 
the shell. We find also an infinite family of "anti-resonances" , where the physical radius of 
the shell is constant (to first order). 

The fact that we have an infinite set of modes suggests that these modes could be used 
to solve the initial value problem (in the linearized approximation). In fact, we can formally 
write an arbitrary solution of the field equations as an integral transform involving Bessel 
functions. Unfortunately, we have not found a way to invert this transforms in such a 
way that they can be written in terms of arbitrary initial data, although this seems to be 
possible in principle. We discussed the reasons that make this problem special, and provided 
a particular example to illustrate some features of the general solution. 
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